Phase-Controlled Force and Magnetization Oscillations 
in Superconducting Ballistic Nanowires. 
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The emergence of superconductivity-induced phase- 
controlled forces in the (10~^ — 10~^) nN range, and of mag- 
netization oscillations, in nanowire junctions, is discussed. A 
giant magnetic response to applied weak magnetic fields, is 
predicted in the ballistic Josephson junction formed by a su- 
perconducting tip and a surface, bridged by a normal metal 
nanowire where Andreev states form. 
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Formation, energetics and mechanical properties of 
metallic nanowires (NWs) , were predicted in early simu- 
lations [1], and they have been the subject of subsequent 
significant research endeavors [2] . For normal metals the 
oscillatory behavior of the elongation forces [1] with the 
size of a NW formed between a surface and a retract- 
ing tip have been shown to be correlated with a quan- 
tized staircase behavior of the electrical conductance [3- 
7]. However, the influence of superconductivity (SC) on 
the nanomechanical properties of such NWs has not been 
explored yet, and these effects are the focus of this Letter. 

In normal metals the cohesive force in an atomic-scale 
metallic contact can be estimated as Fn Sf/^f , where 
Ef and Xf are the Fermi energy and wavelength. The 
onset of SC introduces a new energy scale, i.e. the su- 
perconducting gap A <C ep, and a new lengthscale, i.e. 
the SC coherence length = T^wf/ttA . On first sight, 
the resulting SC- induced forces are expected to be of the 
order of Fgc ~ A/^Oi and when added to the aforemen- 
tioned normal-metal forces (F„ , which are of the order of 
several nN ), they are estimated to be below the atomic 
force microscopy (AFM) detection limit [8]. However, 
in the superconducting regime, under certain conditions 
all the transverse channels {N±) supported by the junc- 
tion will contribute coherently to the free energy, and 
when N± ^ 1 the above consideration may result in 
a gross underestimation of F^c- In particular, we pre- 
dict that under favorable conditions, in a nanowire con- 
necting two superconducting electrodes, modelled here as 
a transparent short [L <^ ^o) superconducting-normal- 
superconducting (SNS) junction, a (measurable) force 
Fsc ^ (£/^o)(£f/A_f) would be manifested - we refer to 
these forces as Andreev forces (AF); in long SNS junc- 
tions the AF are expected to be significantly smaller. 

To calculate the SC - induced contribution to the total 



elongation force we model the junction comprised of the 
superconducting tip and surface and of the NW bridging 
them, as an SNS junction (the NW remains normal due to 
the suppression of superconductivity in very small struc- 
tures [9]). Andreev reflections occur at the SN bound- 
aries of the junction, and the dimensions of the NW are 
such that transport through it is ballistic. This holds 
true for short SNS junctions, while in long ones impurity 
scattering may occur, resulting in a reduced transmission 
probability. 

The force equals the spatial derivative of the grand- 
canonical potential ^l, i.e. F = —dil/dL. In super- 
conducting junctions there are two contributions to the 
AF, Fa = F^ + Fl ; F±_ is related to the dependence of 
the number of quantized transverse channels, N±, on the 
degree of elongation, and Fj^ originates from the depen- 
dence of the Andreev bound states on the length of the 
junction. 

Short ( L <C ^oj junctions. In the following we fo- 
cus on the SC-induced force oscillations in short junc- 
tions. Since here the Andreev states are independent 
of the mode index [10], the potential can be written as 
ris(iy9, L) w iVj_(I/)n^((^, L), where VIa{v^ L) corresponds 
to a single-channel SNS junction. In a cylindrical geom- 
etry N±{L) ~ ttV/X^L, where V is the volume ( which 
remains constant during the elongation process [1] ). 

In general, both bound (superscript (b) below) and 
scattering Andreev states contribute to f2yi(</', L); for our 
purpose only the bound states are important (see below). 
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The spectrum of Andreev bound states in a single-mode 
ballistic SNS junction of length L and transparency D < 
1 is found from the equation (see e.g. Ref.[ll]) 
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where A^ = Tlvf/L^ D + R — 1, and (p is the su- 
perconducting phase difference. The distinction between 
different channels enters through A^, and it is neglected 
for short junctions. A reduced NW transparency reflects 
the presence of interfacial barriers, such as those likely 
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to form at tip-surface contacts. From Eq.(2) we obtain 
in the limit L <C the bound-state energies 



(3) 



EW = -4-^.AH^M(^l-^^|sin||) , 

where Wif) = yJl-D sin^ f . Using Eq.(l), with CIa 
taken simply as the thermodynamic potential for a two- 
level {Eq'^\ Eq~^) system, yields 
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For low transparency junctions {D <C 1) F^_^^ may be 
approximated as 
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and since F^^' ~ N^_{A/£,o)\/D it is evident that | > 
provided that D ^ (L/^o)^. In contrast, in low 



transparency jimctions D <C (L/^o)^ ^ 1 and the 
contribution dominates. 

For point contacts (i.e. extremely short junctions) 
the AF can be calculated (when D ^ (L/^o)^ ) by tak- 
ing the limit L/^o for the Andreev bound states 
(Eq.(3)). In this case, continuum states do not affect the 
free energy (see e.g. Rcf.[ll]) and, therefore, they do not 
contribute to the force. In contrast, for low transparency 
junctions the i/^o-corrections to the bound state ener- 
gies determine the force oscillations. To this order the 
continuum states do contribute to the free energy and 
they can change the dependence of the AF on the phase 
difference. Thus, Eq.(5) can be considered as an esti- 
mate of the AF in a short junction, the contribution 
of the continuum states (F^^) is extremely small, i.e. 

F['\D<^l)^{L/i;of{eF/XF). 

The phase-dependent force in a superconducting quan- 
tum point contact (QPC) (I? = 1) is related to the quan- 
tized Josephson current Jg (see Ref.[10]) 
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It is the latter, phase-dependent, term that is directly 
related to the quantized Josephson current. 

The amplitude of the AF oscillations may be readily 
estimated as follows: for D ^ 1 and L the am- 

plitude of the Andreev force is of the order of F^^^ ~ 
N±A/L ~ Lef/^oXf ~ {L/(,o)[nN]; in the ballistic 
regime for a non-transition metal ^ 10^"'' — IQ^^cm. 
Using state-of-the-art instrumentation such forces (e.g. 
10"^ - lQ~^nN), can be measured [8]. 

Long {L 3> ^o) junctions. For a long transparent 
{D = 1) junction the spectrum of the Andreev-Kulik 
(AK) levels \En\ < A takes the form [12] 

i;±~7rAi(^±£ + ^+n), n = 0,±l,±2,... (8) 

Since this spectrum does not depend on the supercon- 
ducting gap A, for temperatures T <C A all the thermo- 
dynamic properties of a long ballistic junction are essen- 
tially material independent. Evaluation of the Josephson 
current in this case is equivalent to the calculation of the 
persistent current for chiral fermions on a ring with cir- 
cumference 2L [13] . The corresponding phase-dependent 
part of the thermodynamic potential f2^((^) can be eval- 
uated for the AK spectrum, yielding 
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The force oscillations induced by the AK level struc- 
ture in a single-channel long SNS junction are (|<^| < tt) 



At 
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Fa ^ { (10) 
.167r-^ exp (- 1^) cos</p , T > Al . 

We note that the above is equivalent to the Casimir force 
[14]. We focus here only on the phase-dependent part of 
the thermodynamic potential, fJ^ (</?), and the resulting 
Andreev (or Casimir) force Fa . since, as aforementioned, 
the force in superconducting junctions (Eq.(lO)) is added 
to a much larger phase- independent term (~ £f/Af) that 
dominates the cohesive force in metallic NWs. 

In a multi-channel junction the thermodynamic poten- 
tial is the sum over transverse channels {In) 



(11) 



The force oscillations (portrayed by the dependence of 
the force on the contact size) are determined by two dis- 
tinct contributions: (i) a large phase-independent term 
( operative also in normal-metal NWs ) of the order of 
N±Sf/Xf originating from incoherent contributions of all 
the conducting electrons to the thermodynamic potential 
[4-6], and (ii) a coherent SC- induced force (Eq.(4)-(6)). 



where i}^^\p) is given by Eq.(9) with A^ — ,oup 
substituted for A^. For a long junction the Fermi veloc- 
ity enters explicitly the expression for a single channel 
supercurrent, and the total current in a multi-channel 
junction strongly depends on the junction geometry [15]. 

We will model the normal part of a long SNS junction 
by a cylinder of length L and cross-section area S = V/L. 
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Assuming hard-wall boundary conditions, the electron 
longitudinal velocity in the (In)— th channel is 



Ma 



[L) = ^2{ep-a^fJ/m, 
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where 7;„ are the Bessel function zeroes: J; (7i„) = 0, 
and a = h'^TrL/2mV. The dependencies of the AF on the 
phase difference {(p) and on the length (L) of the NW are 
displayed, respectively, in Figs.l and 2, where we show 
AF{(p) = FAif) — i^/i(0). From Fig.l we observe that 
the force is enhanced at special values of the phase dif- 
ference (fir = 7r(2r + 1), (r = 0, ±1, ±2, ...). At tp = (pr 
one of the AK bound states coincides with the Fermi 
energy and, most significantly, this state is 4iVj_-fold de- 
generate [16], thus amplifying its contribution. Direct 
observation of the SC-induccd nanomechanical effect pre- 
dicted here may be obtained through : (i) generation of 
a NW of length L via separation of an AFM tip-surface 
contact, using a superconducting material (e.g. Pb) at 
T < Tc, followed by (ii) measurement of the force re- 
quired to maintain the NW length (L) as a function of 
variations of the phase-difference across the SNS junction 
(as seen from Fig.l this force maximizes at = tt ). 

The variation of the elongation force (for (p = w) with 
the NW length is shown in Fig. 2. We note first that even 
though the number of open channels is very large for the 
NW junction shown in Fig. 2, the magnitude of the forces 
is significantly smaller than in the case of short junctions 
(see previous subsection) [17]. The aperiodic variations of 
the AF originating from the change in the number of open 
channels upon elongation, are particularly pronounced at 
lower temperatures. Note however, that such aperiodic 
variations occur also for normal metal NW [4,6] and con- 
sequently separation of the SC-induced contribution may 
be difficult. 

Next we consider a multichannel superconducting long 

junction in a weak magnetic field iibB <^ A (where /is 
is the Bohr magneton), apphed locally (i.e. only to the 
normal metal nanowire part of the SNS junction). Here, 
the only [18] influence on the AK levels is through the 
Zeeman coupling of the electron spin s to the magnetic 
field, Hz = 5/iBS-B {g is the g- factor). The thermody- 
namic potential SilA{p,B) = ^Ia{^,B) — il^(0,0) takes 
the form (see. Eq.(9)) 
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where x = ^z/^l and Az = qi^bB is the Zeeman 
energy splitting. Note that the infiuence of the Zeeman 
interaction on the thermodynamics of the SNS junction 
is equivalent to the influence of a gate voltage on the 
thermodynamic properties of quantum rings [19]. 

The magnetization Ma = —d5rtA{p,B)/dB at low 
(T <C Al) and high (T > A^) temperatures is given 
for a single channel junction as 



-giiB (f ) , T = ; Ixl < TT 
[SgnB-^e cos ((£3 sin (x) , T > Al . 
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Note, that the SC-induced magnetization Ma, can be 
of the order of several ^b (if 5 ^ 1) even for a single- 
channel junction, and it is insensitive to the supercon- 
ducting phase difference at low temperatures. For a mul- 
tichannel quantum junction at low temperatures the de- 
pendence of Ma{p) exhibits typical resonant behavior 
at the resonant phases pr, as shown in Fig. 3. This is a 
manifestation of the effect of "giant oscillations", known 
previously for conductance oscillations [18]. At these 
phases Andreev states of energies Ea = ±giJ.BsB become 
2A/'j_-fold degenerate [20], leading to giant enhancement 
of thermodynamic and kinetic characteristics of ballistic 
junctions in magnetic fields. 

Since at resonance the coherent contribution (oc A''j_) 
of all transverse modes dominates the magnetization, we 
predict at low temperatures {T <C Az): (i) a giant re- 
sponse (oc N±) of an SNS junction to a magnetic field, 
and (ii) a step-like behaviour of the magnetization as a 
function of the wire diameter. At other values of the 
phase difference, different transverse channels contribute 
to SQa with different periods (i.e. in general, incoher- 
ently), resulting in a complex structure of the magnetic 
oscillations. 

In most cases a supercurrent is suppressed by the 
Zeeman interaction [21]. A magnetic field would also 
suppress the predicted Andreev force SFa(p,B) = 
—dSilAip, B)/dL. At low temperatures (T ^ Al) the 
force (which is periodic both in the phase, p, and in the 
dimensionless Zeeman energy splitting x = Az / A^) can 
be written for a single-channel junction as ( |<^|, |x| < ti" ): 
SFa ^ {AL/2nL) [{^)^ - (x)'] • 

In summary, we predicted and illustrated that super- 
conductivity induces in quantum wires phase-dependent 
forces correlated with the supercurrent. At resonance 
values of the superconducting phase difference these An- 
dreev forces become measurable; (nN scale;). Further- 
more, we predict giant magnetization (of the order of 
N±fiB) of ballistic SNS junctions in a weak magnetic field 
at low temperatures T <C Az- Since low-temperature 
STM with superconducting tips has been already demon- 
strated [22, 23] and used to form Josephson junction 
[24,25], and in light of improved force-detection capa- 
bilities (extending to 10~^ — 10~^nA'')[8], the above pre- 
dictions provide the impetus for future experiments. 
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Fig. 1 Krive et. al. 

FIG. 1. AF{ip) = Fa{<p) - Fa{0) vs the phase difference 
for a NW SNS junction. We use the physical parameters 
for Pb, i.e. the Fermi energy is ef = 1.5 x 10~^^erg, and 
VF = 1.83 X 10*cm/s. The volume = 5 x lO"^'^ cm^, and 
the length of the junction L = 10~'*cm. Results are shown 
for two temperatures, both below Tc{Pb) = 7.18K. The force 
was calculated as follows: first, using Eq.(9), the grand canon- 
ical potential was found for each transverse mode (with dif- 
ferent values of v''p and, therefore, different = hv^p I 
see Eq.(12)). The total potential of the junction is the sum 
of contributions of all the transverse modes (Eq.(ll)). The 
derivative of the potential with respect to the length of the 
junction was evaluated numerically. 
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Fig. 2 Krive et. ai. 

FIG. 2. AF((y3 = 7r) vs the length of the junction, cal- 
culated for different temperatures. The parameters of the 
junction and the method of calculation are as in Fig.l. 
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Fig. 3 Krive et. al. 
FIG. 3. Magnetization of the junction in a magnetic field 
of 5 Oe, plotted vs the phase difference (j) for several tem- 
peratures. The parameters of the junction are as in Fig.l. 
The magnetization was evaluated as a numerical derivative of 
dQ{(p,B) (see Eq.(13)) with respect to B. 
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